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The finite-T DMRG method is applied to the one-dimensional Kondo lattice model to calculate 
dynamic correlation functions. Dynamic spin and charge correlations, St(cj), S c (u}), and N c (u>), and 
quasiparticle density of states p(ui) are calculated in the paramagnetic metallic phase for various 
temperatures and hole densities. Near half filling, it is shown that a pseudogap grows in these 
dynamic correlation functions below the crossover temperature characterized by the spin gap at half 
filling. A sharp peak at to — evolves at low temperatures in S{(u>) and N c (ui). This may be an 
evidence of the formation of the collective excitations, and this confirms that the metallic phase is 
a Tomonaga-Luttinger liquid in the low temperature limit. 



Heavy fermion systems have attracted much attention 
for more than a decade because of their enormous mass 
renormalization and diverse ground states including un- 
conventional superconductivity. Q The periodic Ander- 
son model (PAM) and the Kondo lattice model (KLM) 
are their canonical models, and their properties have in- 
tensively been investigated. Q In the framework of the 
PAM, hybridized band picture with strong renormaliza- 
tion is widely accepted for a scenario of the formation 
of heavy quasiparticles. 3j|(| The formation of heavy 
quasiparticles is supported by numerical studies on the 
PAM in infinite dimensions (D=oo) |Qj8) and in one di- 
mension (ID), [jj 10 1 On the other hand it is a nontrivial 
question how the heavy fermion state is formed in the 
KLM, |ll]] since the charge degrees of freedom are com- 
pletely suppressed for the localized /-electrons. For this 
problem it is important to understand the relation be- 
tween the Kondo effect and the hybridization picture. 
An important energy scale is the coherence tempera- 
ture of the heavy quasiparticles, T co h, but several ap- 
proaches seem to give controversial results. Gutz wilier 
approximation suggests that T co h is enhanced compared 
with the single-impurity Kondo temperature Tk, 0] while 
Nozieres predicted that T co h is rather suppressed due to 
the exhaustion mechanism, [fl2|| Several numerical calcu- 
lations on the PAM in D=oo and ID suggest the sup- 
pression of T coh . 

Intersite correlations are believed to play an important 
role for the formation of heavy quasiparticles in real com- 
pounds, but they are not taken into account in D=oo. ID 
systems provide complementary information, since the 
intersite correlations generally become more dominant in 
lower dimensions. Precise numerical calculations are also 
feasible in ID compared with 2D or 3D. Therefore the 
ID KLM may be another appropriate starting point to 
obtain reliable information on dynamics, although the 
ground state may not be a Fermi liquid. 

Ground-state properties of the ID KLM have been 
studied extensively and its magnetic phase diagram is 
now determined. Jll|[l4[ At half filling the ground state is 
a Kondo spin- liquid insulator with a spin gap A s , 1 15 2(J 
and the charge gap A c is much larger than A s due to the 



correlation effects on the gap formation. [[l7]|lj| Upon fi- 
nite hole doping 5, both spin and charge gaps close and 
the ground state is expected to belong to the universality 
class of the Tomonaga-Luttinger liquid (TLL). pl|-p3[ 

Thermodynamic properties at finite S have been stud- 
ied by the finite-temperature density-matrix 
renormalization-group (finite-T DMRG) method (2J |||] 
and it is found that there exist two crossover tempera- 
tures. |2^] The first crossover is observed as a broad peak 
in the T-dependence of the spin susceptibility x s when 
the temperature is lowered from infinity. The crossover 
temperature Tj* may be defined by its peak position, and 
calculations for various J's show that Tj* is scaled by A s . 
p9| , ^8| The charge susceptibility Xc and Xs both decrease 
below T* but they turn to increasing again at further low 
temperatures. The susceptibilities arc expected to finally 
approach the T = value determined by the velocities of 
the collective excitations of the TLL ground state, and 
this saturation is actually observed in Xc for some J and 
S. [^8| The second crossover temperature T 2 * may be de- 
fined to characterize this saturation, and is essentially 
determined by either spin or charge velocity depending 
on the channel. The zero temperature susceptibilities cal- 
culated by the DMRG method show a diverging behavior 
as S — > 0, and this means vanishing of T 2 *. |28| 

In the present work we calculate temperature and dop- 
ing dependence of the dynamic correlation functions and 
clarify the character of these crossovers. It will be shown 
that the higher crossover temperature T* corresponds to 
the characteristic temperature of pseudogap formation 
in the dynamic spin and charge structure factors. At the 
same time a sharp peak structure appears at uj = 0, in- 
dicating the formation of the collective excitations of the 
TLL at low temperatures. The lower crossover temper- 
ature T 2 may correspond to the coherence temperature 
of these collective excitations. 

The Hamiltonian of the ID KLM is described as 



* = -*£ 



cT a +u + H. 



i,s,s' 



' c is c is' 



(i) 
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with standard notations. The density of conduction elec- 
tron n c is unity at half filling, and hole doping (n c = 1—5) 
is physically equivalent to electron doping (n c = 1 + 6) 
due to particle-hole symmetry. 
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FIG. 1. Dynamic spin correlations of the /-spins in the 
ID KLM. A s = OAt is the spin gap at 5 = 0. 
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FIG. 2. Dynamic spin correlations of the conduction elec- 
trons. 
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FIG. 3. Dynamic charge correlations. A c = 1.4t is the 
charge gap at 5 = 0. 

Dynamic correlation functions of local quantities can 
be calculated at finite temperatures for infinite-size sys- 



tems by the DMRG method |23] applied to the quantum 
transfer matrix. []30f Imaginary-time correlation func- 
tions are first calculated from the eigenvector of the max- 
imum eigenvalue of the quantum transfer matrix, |3l|-|35|l 
and then they are analytically continued to real frequency 
using the maximum entropy method. The advantages 
of this method are that the finite-T DMRG method has 
no statistical errors and that we do not need the extrapo- 
lation on the system size. This approach was first applied 
to the insulating phase of the ID KLM and the many- 
body nature of the gap formation is revealed. |3l],[52| The 
present study is the first application to a metallic phase. 
In our calculations we usually keep 50 states in the finitc- 
T DMRG procedure with the Trotter number 60. 

We have calculated temperature dependence of several 
dynamic correlation functions. The results of the local 
spin dynamics of the /-spins, Sf(u>) = J ^Sf(q, lu), are 
shown for J/t = 1.6 and 5 = 0.2 in the inset of Fig. [I]. 
Note that the complete suppression of charge fluctuations 
for /-spins imposes the sum rule, / duS{(ui) = 1/4, inde- 
pendent of J, 5 and T. We can see that a peak structure 
appears around lu <~ A s = 0.4t at T < 0.2t, and similar 
peak structure is also observed for different Kondo cou- 
pling J/t = 1.2 at T < 0.06t, where A s = 0.16*. Based 
on these results, we may conclude that the characteristic 
temperature of the peak formation is scaled by the higher 
crossover temperature Tj" determined by Xs(T). At the 
same time another peak structure grows at u> = 0, when 
S is finite. It suggests the formation of the collective spin 
excitations of the TLL at low temperatures. We expect 
that the low energy part finally approaches \m | Kc , which 
is predicted by the TLL theory at T = 0. [|34| 

The doping dependence of S{(lo) is shown at the low 
temperature T = O.OAt <§; A s in the main panel of Fig. [|. 
With increasing 5, the peak intensity at lu = grows, 
while the intensity around uj = A s is reduced as a conse- 
quence of the sum rule discussed before. We have checked 
that the peak intensity at u) = is 0.25(5 within a few per- 
cent. The energy scale estimated from the peak width at 
lu = at T = 0.04i is smaller than the lowest tempera- 
ture in our calculations. This means the weight around 
uj = yields nearly free spin degrees of freedom with den- 
sity 0.255 down to around this temperature T ~ 0.04t. 
This is consistent with the behavior of the static spin sus- 
ceptibility, Xs = ( 5/(4T), observed at low temperatures. 

In contrast to S{(uj) the lu = peak is small 
in the conduction-electron spin correlations, S c (lo) = 
I ?~Sc(q> w), as shown in Fig. ||. This shows that the low- 
energy part of the spin degrees of freedom of the conduc- 
tion electrons are nearly exhausted to form spin singlet 
with /-spins and thus a clear peak is formed at lu = A s . 
Therefore, the low energy spin dynamics is mostly dom- 
inated by the /-spin degrees of freedom. However, the 
intensity of the peak around lu ~ A s is less than 1/3 of the 
corresponding peak in Sf(cj) and over a half of the total 
weight extends over higher frequencies t < lu < 5i. This 
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means that although the low energy part is dynamically 
coupled with /-spins, the majority part of the conduc- 
tion spin degrees of freedom have another energy scale 
of almost the band width ~ At s> A s when the Kondo 
coupling is small J < At. This shows that only the con- 
duction electrons close to the Fermi level screen the /- 
spins as pointed out by Nozieres. jl2j Rather surprisingly, 
even though the screening by conduction electrons is not 
complete, the intensity at u> = in Sf(u>) is almost <5/4, 
just like the J = oo case, where each conduction electron 
screens one /-spin. This shows the importance of the /-/ 
spin correlations for the formation of the Kondo singlet 
state. 

The doping dependence of the dynamic charge corre- 
lations, N c (w) = J ^N c (q,cj), is shown in Fig. |. At 
half filling, the charge excitations are exponentially sup- 
pressed below the crossover temperature T* at u> < A c . 
p9| Upon hole doping, a sharp peak appears at u> = 
as in Sf(uj). This indicates the formation of collective 
charge excitations of the TLL as discussed for S{(ui). The 
peak intensity increases with 6, and this means that the 
effective carrier density of the TLL is strongly renormal- 
ized from n c = 1 — S. The renormalization of carrier is 
naturally explained in the limit of strong J, where each 
conduction electron forms a local singlet with the /-spin 
on the same site. In this limit, effective carriers are in- 
troduced by hole doping, and their main component is 
the unscreened /-spins with density 5, j21j as discussed 
for St(u>). The TLL theory predicts that the low-energy 
asymptotic form of N c (u>) is \uj\ min ( K o' 4K o-i) near w — fj, 
and we expect that the u — peak finally approaches 
this form in the low temperature limit. 

We now consider crossover behavior of the quasiparti- 
cle density of states, p(u>). Figure || shows the temper- 
ature dependence at S = 0.2 and J/t — 1.6. We can 
see that a pseudogap develops just above the Fermi level 
lu = fi below T/t ~ 0.4. Similar behavior is observed also 
for J/t — 1.2 below T/t ~ 0.15. Based on these results, 
we may conclude that the characteristic temperature of 
pseudogap formation is scaled by T* defined from Xs(T). 
This conclusion is consistent with the results obtained at 
half filling. HI! 

The crossover behavior around T-j* may be explained as 
follows. Below Tj* thermal fluctuations of the /-spins are 
substantially suppressed, since the temperature is lower 
than the characteristic energy scale of the /-spin excita- 
tions A s . When 6 is small, the characteristic time scale 
of the dominant part of the /-spins is given by A^T 1 , and 
this is much longer than the time scale of quasiparti- 
cle propagation r qp , since r qp may be determined by the 
inverse of bare hopping energy and charge gap, t^ 1 and 
AjT 1 . Therefore, concerning the quasiparticle excitations, 
the /-/ spin correlations may be assumed to be static. 
Because of their staggered nature in space, these almost 
static /-/ spin correlations induce the opening of a gap 
in p(uj). Of course, as shown in S(u>) and N(u), both 
the /-spins and the conduction electrons have a slow dy- 
namics when S is finite, and therefore the gap discussed 



here is not a real gap but rather a pseudogap for finite S. 
Thus the pseudogap develops below T*. 




FIG. 4. Temperature dependence of the quasiparticle den- 
sity of states. 
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FIG. 5. Doping dependence of the quasiparticle density of 
states. Quasiparticle gap at S = is A qp = 0.7i. 

We next consider doping dependence to study small 
structures in p(to) near the Fermi level uj = p. The re- 
sults for J/t = 1.6 at T/t = 0.04 are shown in Fig. | at 
various dopings. The pseudogap becomes more promi- 
nent with approaching half filling, where the clear quasi- 
particle gap, Aqp/i = 0.7, exists. The sharp peak struc- 
ture just below u> = fx also grows with decreasing S and 
seems to continuously connect to the gap edge structure 
at 5 = 0. 

The nature of the structure near the Fermi level is not 
fully understood yet. One possible scenario is the mean- 
field type argument assuming the /-spin helical SDW 
order with wave number 2k p — 7r(l — S). Band mixing 
induces gap opening at k = tt(±1 + <5)/2, and two new 
van-Hove (vH) divergent singularities appear in each of 
the two hybridized bands. The Fermi level sits between 
these two new singularities in the lower hybridized band, 
as far as 5 is small. A slightly different scenario is that the 
gap is induced by the short-range /-spin correlations with 
wave number n. Then, there appears only one vH singu- 
larity in each band. The peak-like behavior observed in 
Fig. H may be identified as the lower vH singularity in 
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FIG. 6. Doping dependence of the quasiparticle density of 
states. Quasiparticle gap at S — is A qp = 0.47t. 

the first scenario, while in the second one it is a conse- 
quence of the shift of the Fermi level towards the top of 
the lower band as S — > 0. However, it is not straight- 
forward to describe the behavior of p(u>) just above the 
Fermi level in terms of these pictures. Another quite 
different scenario is that this peak is due to the Kondo 
singlet formation and reminiscence of the Kondo reso- 
nance. This may also be expressed as the renormal- 
ized hybridized bands as proposed for the PAM, with 
the Fermi level inside the lower band. The renormal- 
ized hybridized bands are actually observed for the PAM 
in D=oo H and ID. (To) Although the results of the 
ID PAM indicate a quite symmetric p(u>) near u) = fi 
for conduction electrons, which differs from our results, 
there is a considerable asymmetry in the results for the 
D=oo PAM, similar to our results. With decreasing 8, 
the Fermi level shifts towards the top of the lower band, 
and p(u>) near the Fermi level grows accordingly, which 
is consistent with the behavior in Fig. ||. Therefore, this 
scenario is quite promising. Anyway, the slow dynamics 
of /-spin correlations may be important to understand 
the low energy dynamics of quasiparticle motion, and we 
will make further investigations in future study. 

As the temperature decreases down to zero, the in- 
teraction between the effective carriers becomes relevant 
and the low energ y e xcitations are expected to be de- 
scribed as a TLL. [^||,^3) The second crossover tempera- 
ture is defined to characterize this behavior, as discussed 
before. Figure ^ shows p(u>) for the smaller coupling 
J/t = 1.2. In this case, the second crossover temperature 
is relatively high compared with the case of J/t = 1.6, 
p8f and the asymptotic TLL behavior may be observed. 
For 8 = 0.2 and 0.3, a dip structure is now observed 
at the Fermi level to = p, which is absent in Fig. |^. It 
is known that the TLL has p(ui) with a dip structure 
at the Fermi level as p(u>) ~ \oj — p\ a at T = where 
a = (K c - 1) 2 /(4K C ). Q Here K c is the Luttinger liq- 
uid parameter and less than 1/2 in the ID KLM. ]23| ] 
Thus the dip structure in Fig. ^ is consistent with the 
TLL picture. Such dip structure is not found for small 
doping 8 = 0.1. This agrees with the previous study of 



the thermodynamics in the ID KLM, [^8| which shows 
that T 2 * i s lower for smaller 8 and vanishes as 8 — > 0. 

To summarize we have calculated dynamic quantities 
at various temperatures and hole densities, and clari- 
fied characters of the two crossovers in the paramag- 
netic metallic phase. Below the first crossover temper- 
ature T-j* it has been shown that a pseudogap develops in 
the density of states and dynamic correlation functions, 
and S{(lo) and S c (lo) both show a peak structure around 
ui — A s as in the half-filling case. At the same time a 
peak structure appears at lu — in S[(w) and N c (u), and 
its intensity increases with hole density 8. The u) = 
peak indicates that as a consequence of local Kondo sin- 
glet formation effective carriers are strongly renormalized 
to have density 8 and small energy scale. The increase 
of the peak intensity with 8 is naturally explained in the 
limit of strong J, where effective carriers of the TLL are 
the unscreened /-spins whose density is 8. We note that 
this is consistent with the large Fermi surface. J2^j2^| Be- 
low the second crossover temperature T 2 * , the interaction 
between the effective carriers becomes relevant and the 
renormalized carriers are expected to evolve into a TLL. 
This is supported by several sets of the present results. 
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